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Abstract 

We examine the decay of passive scalars with small, but 
non zero, diffusivity in bounded 2D domains. The ve- 
locity fields responsible for advection are smooth (i.e., 
they have bounded gradients) and of a single large scale. 
Moreover, the scale of the velocity field is taken to be 
similar to the size of the entire domain. The importance 
of the initial scale of variation of the scalar field with 
respect to that of the velocity field is strongly empha- 
sized. If these scales are comparable and the velocity 
field is time periodic, we see the formation of a periodic 
scalar eigenmode. The eigenmode is numerically realized 
by means of a deterministic 2D map on a lattice. An- 
alytical justification for the eigenmode is available from 
theorems in the dynamo literature. Weakening the no- 
tion of an eigenmode to mean statistical stationarity, we 
provide numerical evidence that the eigenmode solution 
also holds for aperiodic flows (represented by random 
maps). Turning to the evolution of an initially small 
scale scalar field, we demonstrate the transition from an 
evolving (i.e., non self similar) pdf to a stationary (self 
similar) pdf as the scale of variation of the scalar field 
progresses from being small to being comparable to that 
of the velocity field (and of the domain). Furthermore, 
the non self similar regime itself consists of two stages. 
Both the stages are examined and the coupling between 
diffusion and the distribution of the Finite Time Lya- 
punov Exponents is shown to be responsible for the pdf 
evolution. 

PACS number 47.52. +j 
1. Introduction 

Starting with the work of Batchelor 1 the study of passive 
scalars in smooth velocity fields has been the subject of 
numerous investigations. Originally posed in 3D 1 , the 
problem considered a situation where the viscosity (v) 
of the flow is much greater than the diffusivity (k) of 
the scalar. Therefore, even at large Reynolds numbers, 
for length scales below the viscous cutoff (l u ) and above 
the diffusive cutoff (l K ), the so called Batchelor regime, 
one has a smooth velocity advecting a diffusive tracer. 
Since then the problem has grown to encompass a vari- 
ety of phenomena such as chaotic advection (see for eg. 



Ottino 2 ) and scalars in inverse cascading 2D turbulent 
flows (see Section III in Falkovich et. al. 3 for a recent 
review). Of particular interest to us are applications in 
the realm of geophysical fluid dynamics, specifically, the 
mixing of scalars along isentropic surfaces via large scale 
atmospheric flows 4 , 5 . 

The equation governing the passive scalar ((f) is, 

^ + (u-V)0- K V 2 (1) 

Let us denote the size of the domain by L. This is a 
linear equation for </>, the velocity field is part of the 
prescribed data. The usual conditions imposed on the 
velocity field are that it be divergence free and smooth, 
i.e., V • u = and Vw| < oo everywhere in the domain. 
The questions asked are usually of the form : given an 
initial scalar field is it possible to determine the long 
time behaviour of the solution ? Does one see the emer- 
gence of an eigenmode with a well defined decay rate, in 
other words what is the nature of the spectrum of the 
advection diffusion operator in a suitably defined func- 
tion space ? Usually, given the aperiodic nature of the 
advecting velocity fields the question is better posed in a 
statistical sense, i.e., is the asymptotic pdf of the scalar 
field self similar or is it an ever evolving (non self similar) 
entity ? As things stand, the results found in the liter- 
ature are fairly divided. One of the aims of this paper 
is to unify these results by bringing out certain salient 
features (such as scale separation) which have not been 
emphasized previously. 

Most theoretical studies have focused on the case when 
the initial scale of variation of the scalar field (say l s , 
l s >> l K ) is much smaller than that of the velocity field 
(say l v , we take l v ~ L). The problem then is to deter- 
mine the statistical properties of the scalar field at scales 
between l v and l K . A successful approach in this direc- 
tion has been to shift to a comoving reference frame, use 
the decomposition u a = (d a Ufj)(t) rp = o~ a p(t) rp and 
deal with the the effective equation that results from a 
substitution in Eq. (1) (see for eg. 6 , 7 and the references 
therein). The expectation 7 is that the scalar field will 
have non self similar pdf's. Or, defining the moments 
of the scalar field as < \<fi(x,y,t)\ n >, Balkovsky and 
Fouxon 7 explicitly show that, 

< \4>(x,y,t)\ n >~ e -«™* ; t>T (2) 



and that a n is a nonlinear function of n (T is a diffusive 
time scale which will be clarified later) , implying the non 
self similarity of the pdf's. Recent experiments 8 with 
passive scalars injected at point sources in inverse cas- 
cading 2D turbulent flows seem to validate some of these 
predictions. Concrete evidence is available from numeri- 
cal work demonstrating the nonlinearity of a n for scalars 
advected by realistic atmospheric winds 5 . It should be 
noted that these simulations are over a time scale which 
represent the transient scalar behaviour rather than the 
asymptotic large time solution. 

On the other hand, there is a small body of work which 
deals with the case when l K « l s ~ l v ~ L, i.e., the ini- 
tial scale of variation of the scalar field is comparable to 
that of the velocity field, both of which are in turn com- 
parable to the scale of the domain. 1 The observation, 
based on numerical work 10 , 11 , is that the scalar field en- 
ters an eigenmodc (termed a "strange eigenmode" by 
Pierrehumbert 10 due to its spatial complexity). The pdf 
of the scalar becomes self similar (after a suitable nor- 
malization by the variance) or in terms of the moments, 
a n ~ n. Recently experimental evidence for the emer- 
gence of such an eigenmode has been provided in the 
case of time periodic velocity fields 13 . A line of attack 
on this problem has been to try and relate the Finite 
Time Lyapunov Exponents (FTLE's) of the advecting 
flow to the statistical properties of the scalar field (see 12 
and the references therein). Even though this work 12 is 
successful in describing the initial stages of the problem, 
the FTLE's are by definition local entities and appear to 
be unsuitable precisely when l s ~ l v 2 . 

In this study we work with simple prescribed veloc- 
ity fields in 2D that have a single large scale (of the 
kind encountered in studies of chaotic advection). It 
is worth emphasizing that the power spectrum of such 
velocity fields is discrete. This is in contrast to techni- 
cally smooth velocity fields which possess a heirarchy of 
scales. In other words the velocity fields considered are 
Lipschitz continuous over all scales of interest (i.e., be- 
tween l K and l v ~ L). The large single scale, divergence 



1 We always take l v ~ L. There are other problems of in- 
terest where l v « L which are studied in the context of 
turbulent diffusion (see for eg. the comprehensive review by 
Majda and Kramer 9 ). 

2 Note : A recent paper 34 analyzes the work of Antonsen 
et. al with an aim to understand the exponential decay of 
scalar moments. This paper, Wonhas and Vassilicos 34 , iden- 
tifies a "global" mechanism responsible for the exponential 
decay and also hints at the importance of the scale separa- 
tion we have emphasized in this paper. 



free and smooth nature of the velocity field implies that 
the trajectory equations do not show explosive (implo- 
sive) separation (collapse) as is generally expected in tur- 
bulent (or multiple scale) incompressible (compressible) 
velocity fields (see for eg. Gawedzki and Vergassola 14 ). 

Our first aim is to give further justification for the 
eigenmode in periodic and aperiodic velocity fields. We 
explicitly construct an area preserving deterministic 2D 
map which, in spite of having weak barriers, demon- 
strates the emergence of a periodic eigenmode. More- 
over, we point out theorems in the dynamo literature 
which are applicable to the passive scalar problem, 
thereby providing analytical justification (in terms of the 
purely discrete spectrum of the advection diffusion op- 
erator) for the eigenmodes. The aperiodic case is free of 
barriers due to the random nature of the maps employed 
but is slightly novel in the sense that the notion of an 
eigenmode has to be weakened by invoking statistical 
stationarity (i.e. self similarity of the pdf). Secondly, 
we present a unified picture of the two seemingly dis- 
parate cases mentioned above. Starting with a scalar 
field which satisfies l s << l v ~ L, we observe an evo- 
lution in the shape of the pdf (i.e., the pdf is non self 
similar). This evolution ceases when the scales of the 
velocity and scalar fields become comparable (i.e., the 
field enters an eigenmode with a self similar pdf). The 
various regimes that the scalar field encounters are an- 
alyzed and a simple example, motivated by the work of 
Balkovsky and Fouxon 7 , is used to elucidate the cou- 
pling between diffusion and the distribution of FTLE's 
and its role in determining the pdf evolution. A discus- 
sion of the regimes, their universality and the issue of 
scale separation concludes the paper. 



2. The Eigenmode 

A. The Eigenvalue Problem 

Physically when l s ~ l v ~ L the scalar field feels the 
effect of the finite size of the domain, moreover, as l s ~ l v 
it is difficult to justify any linearization of the underlying 
nonlinear trajectory equations. Fortunately, the linear 
nature of Eq. (1) lends itself to an eigenvalue analysis. 
Recently Fereday et. al 15 have looked at the eigenvalue 
problem for a specific velocity field (represented by the 
Bakers map). We consider the eigenvalue problem in its 
full generality as has been done for the magnetic field 
in the kinematic dynamo literature 16 , 17 , 18 . Most of the 
material in this subsection is well known in the dynamo 
literature, we present it here merely to put things in a 
well defined framework. Let us recast Eq. (I) as, 

^=£ K ; U = «V^-(«-V)^ (3) 



For a steady velocity field one can separate time by as- 
suming a solution of the form <j>(x,y,t) — 4>(x, y)e at . 
Hence the eigenvalue problem is, 

= <t0 (4) 

Here <j> e 5(D), where 5(D) is a Banach space of square 
integrable functions 3 defined on the domain D. We take 
D to be [0, 2w] x [0, 2ir] with opposite sides identified (i.e., 
a torus). The object of interest is the spectrum of the 
operator C K acting in 5(D). As we are working in a Ba- 
nach space of infinite dimensions the spectrum depends 
strongly on the nature of the operator 19 . 

If the operator in question were compact 4 then the 
spectral theory of finite dimensional operators would 
carry over to our infinite dimensional problem 19 . Unfor- 
tunately differential operators are usually not compact, 
moreover they are unbounded. One of the techniques 
used to handle such operators is to gain control over 
their resolvents (see Kato 19 Chap. 3), the case when the 
operator is closed and its resolvent is compact is partic- 
ularly clear cut 19 . As it turns out this is precisely what 
happens to C K for steady velocity fields and non zero k, 
i.e., it is (weakly) closed on 5(D) and has a compact 
resolvent 20 . 5 Hence its spectrum consists of a discrete 
set of eigenvalues of finite multiplicity (Kato 19 p. 187), 
in fact it has a complete set of cigenfunctions 20 . Phys- 
ically the implication of a purely discrete spectrum is 
that the eigenvalues are well separated, as time goes on 



3 For finite times, given that |Vw| < oo, we are guaranteed 
that the scalar field will remain square integrable. 

4 Consider an infinite sequence of functions ipi,ip2--- in a 
complete function space. Let us denote the action of an op- 
erator A on these functions by the sequence Aipi , Aifa.... If 
A is closed with respect to this space and further if one can 
extract a convergent subsequence from Aipi, Alfa... then A is 
called a compact operator. 

5 The cited reference 20 deals with a similar operator which 
appears in the magnetic dynamo theory, specifically the op- 
erator they work with is, 

M K B — k\7 2 B — (u ■ V)£? + B ■ Vu (5) 

where B is the magnetic field. The difference in C K and 
M K is the additional "stretching" term in M K - As this term 
does not contain derivatives of B it is quite intuitive that the 
spectral properties of M K should carry over to £ K . In fact, 
it is quite reasonable that as long as the gradients of the 
velocity field are bounded, the spectral properties of both C K 
and M K follow those of the Laplacian as the highest order 
derivatives are contained in the Laplacian term. 



4>(x, y) will assume the form of the eigenfunction corre- 
sponding to the largest eigenvalue. In the more general 
case of a periodic velocity field we have to consider the 
Floquet problem involving the propagation of (j>(x,y,t) 
to 4>{x,y,t + T) where T is the periodicity of the un- 
derlying flow. Denoting the propagator by T K (T), the 
eigenvalue problem is of the form, 

%{T]4, = e? T '4> (6) 

The discreteness of the spectrum of T K persists but 
the existence of a complete set of eigenfunctions is not 
gauranteed 20 . Essentially when the flow is smooth and 
k > 0, we are assuming that the spectral properties of 
C K follow those of A4 K (Eq. (5)) for which there exist 
rigorous results in the cited dynamo literature. So, from 
a mathematical point of view the problem has a purely 
discrete set of eigenvalues and a (possibly incomplete) set 
of eigenfunctions for both steady and time periodic ve- 
locity fields. In passing we mention that the k — > limit 
is expected to be quite delicate; interesting aspects not 
encountered in the k ^ problem are likely to appear 
in this limiting process (see Bayly 17 for some illuminat- 
ing examples in regard to the analogous scalar dynamo 
problem) . 



B. Deterministic Maps (Periodic Velocity Fields) 

We proceed to see if we can numerically realize the eigen- 
functions corresponding to the above mentioned eigen- 
values. Chaotic maps on a lattice are utilized to rep- 
resent the advection process (see Pierrehumbert 11 for a 
succinct overview of lattice maps in advection). The 
advantage of the lattice maps is the exact preservation 
of scalar moments and their numerical efficiency. This 
"pulsed" advection is followed by a diffusive step (for 
eg. see 11 or 10 ). Usually maps with random phase shifts 
at each iteration are employed so as to break the KAM 
tori which are generically present in deterministic area 
preserving maps. As this procedure gives an aperiodic 
time dependence to the flow it is unacceptable in the 
eigenvalue formulation outlined in the previous section. 
Unfortunately, we are unaware of any area preserving, 
nonlinear and continuous map which has been proven to 
be mixing 6 over the whole domain D. The only way we 
have of seeing that there are no regions where the scalar 
field remains trapped is to do so numerically, i.e., to ex- 
amine the variance as a function of time. The particular 
map we use mimics an alternating nonlinear shear flow, 



6 Mixing in the sense of dynamical systems (see for eg. 21 ). 



x n +i = x n + A x sin(Si2/„ + si) + A 2 sm(B 2 y n + si) 
J/n+i = J/n + Ci sin(£>ia: n+ i + s 2 ) + 
C 2 sin(L) 2 a;„ + i + s 2 ) 

Here si, s 2 are constant random numbers used to get the 
sine functions out of phase and x n ,y n are mod(0, 2ir). 
The constants are A\ = 4, B\ = 1, A 2 = 1, B 2 = n, C\ = 
1,-Di = 1,C 2 = j,.D 2 = § and the initial condition of 
the scalar field is <j>(x,y) = cos(x)cos{y) (note that the 
scale of the scalar field is comparable to the scale of the 
flow). The diffusive step is implemented as, 

t/Hj = (I - D)<Pij + A&j (8) 

where Afaj is, 

^4>i,j = J D (<Pi+l,j + </>i-l,j + <Pi,j+l + 4>i,j-l) (9) 

Here < D < 1 is the diffusion coeffecient and i, j are 
the indices of the lattice. As is expected, for certain 
combinations of si, s 2 there are persistent KAM barriers 
and we do not observe a decay of the scalar variance as 
would be expected from a globally mixing system. In 
fact, after a short time the diffusive exchange across the 
barriers controls the variance, si = 0, s 2 = provides an 
example of this phenomenon, Fig. 1 shows the variance 
as a function of time for D = 0.5,0.4,0.3 and D = 0.2, 
the control of the variance by the diffusion coefficient is 
evident, moreover, a certain amount of scalar is trapped 
within the barriers (the trapping would be permanent as 
k — > 0, though the scaling of the amount that diffuses as 
a function of k is expected to be nontrivial 22 ). 

On the other hand we get the desired global decay 
(down to machine precision) accompanied by the emer- 
gence of a periodic eigenmode for other pairs of Si,s 2 . 
(namely, si = 5.9698, s 2 = 1.4523 and D = 0.5). Apart 
from using <p(x, y) = cos(a;)cos(j/) as an initial condition 
we also carried out the simulation using a checkerboard 
type initial condition where the field is discontinous but 
square integrable (the initial conditions can be seen in 
the upper panel of Fig. 2). The spatial structure of 
the eigenmode at iteration 250 and 650 for the first ini- 
tial condition can be seen in the first column of Fig. 2. 
The second column of Fig. 2 shows the evolution of the 
scalar field for the second initial condition (the eigen- 
mode is shown at iteration 250 and iteration 650). The 
upper panel Fig. 3 shows a log plot of the moments as 
a function of iteration. 7 The perfectly normal scaling, 



The higher moments involve a fair degree of numerical un- 
certainly, even though they do appear to behave as expected, 
we feel it is better to rely on the pdf 's themselves. 



i.e. a n ~ n, (refer Eq. (2)) is shown in the middle panel 
of Fig. 3. Finally the self similar pdf's in steps of 25 
iterations from iteration 250 to 650 can be seen in the 
lowermost panel of Fig. 3. Fig. 4 shows the same en- 
tities for the second initial condition (see figure caption 
for details). The bimodal nature of the pdf indicates 
that there still exist regions which remain relatively iso- 
lated from each other. In a sense these regions are "large 
enough" to permit chaotic advection of the scalar within 
themselves and the barriers are "leaky enough" so the 
exchange of the scalar across them is strong enough so 
as not to be a controlling factor. Note that due to the 
existence of these barriers the value of D has to be rea- 
sonably large (in our experiments the results for a n were 
identical for D > 0.2), if the system was mixing without 
any barriers then the results would be valid for infinites- 
imally small diffusivity. 

Though the emphasis in the dynamo literature has 
been on the infinite magnetic Reynolds number limit, 
there have been a few studies on the form of the eigen- 
functions at finite magnetic Reynolds numbers. The 
maps used in these studies are chosen for their math- 
ematical properties (such as the Bakers map in 23 or 
the Cat map on an unusual Reimannian manifold in 24 ), 
i.e., one has good mixing properties over the whole do- 
main. Also, these maps allow analytical estimates of the 
form of the eigenfunctions, hence providing some insight 
into the generically singular infinite magnetic Reynolds 
number limit (in terms of convergence to some eigen- 
distribution 18 ). Unfortunately such maps are slightly 
unphysical, their chaotic properties are due to bound- 
ary conditions or due to the curvature of the underlying 
manifold rather than any inherent nonlincarity in the 
maps themselves. Even though we choose not to use 
such maps for advective purposes, we would like to men- 
tion that their results (especially 23 ) are consistent with 
our work. 



C. Random maps (Aperiodic Velocity Fields) 

In the aperiodic case there is no clean way to separate 
time from space as was done in the earlier situations. 
Due to this we weaken the notion of an eigenfunction 
by invoking statistical stationarity, i.e., the scalar field 
is said to enter an eigenmode when the shape of its pdf 
remains invariant with time (upon renormalization by its 
variance as the overall field strength is decaying). Note 
that the earlier classical eigenmode also has self-similar 
pdf's. Of course it has a stronger form of convergence 
where the scalar field itself approaches a stationary (or 
periodic) spatial structure. The random map used to 
represent the aperiodic flow is again a nonlinear shear 
flow, 



x n +i =x n +4 sm(y n + p n ) 
Vn+i = Vn + sin(x rl+ i + q n ) 



(10) 



Now Pn,q n (£ [0j27r]) are random numbers which are 
selected at the beg inning of each iteration. x n , y n are 
mod(0, 2tt), D = 0.25 and the initial scalar field is 
4>{x, y) — cos(0.5x)cos(0.5?/). The map is iterated till the 
variance goes down to machine precision. The variance 
as a function of time can be seen in the upper panel of 
Fig. 5. After an initial transient (which lasts for about 
10 iterations, sec 12 ) the variance decays exponentially. 
The higher moments can be seen in the lower panel of 
Fig. 5. The non- anomalous nature of the exponents 
can be seen in the upper panel of Fig. 6, the implied 
self-similar behaviour of the pdf of the scalar field is 
demonstrated in the lower panel of Fig. 6. The pdf is 
unimodal as there are no isolated regions in the domain 
(the randomness at each iteration is responsible for the 
destruction of any barriers that might exist in the steady 
map). Moreover, the shape of the pdf is characterized 
by a Gaussian core and stretched exponential tails, as 
per the higher resolution studies in Picrrehumbcrt 11 . 



As (f) S B(D), via Fourier's theorem, we can represent 



as, 



cj>(r,t) 



/ 4>( k o, 

Jko 



t) e^^dko 



(11) 



Substituting a plane wave solution in Eq. (1) and equat- 
ing the real and imaginary parts, we get, 



where k{t) is given by, 
- dk(t) 



dt 



= -u- k(t); k(0) = k 



(12) 



(13) 



A. Small time scales 



D. Physical Interpretation 

The physical picture that goes along with the eigenmodc 
is as follows : Consider a filament of scalar whose length 
is comparable to the scale of variation of the velocity field 
(or the typical size of an eddy in the velocity field) . Due 
to the assumed chaotic properties of a generic time de- 
pendent 2D flow, the filament will tend to be stretched 
out. The point is that, as the scale of the filament is 
already the same as that of the flow, instead of being 
merely stretched the filament will fold and start to fill the 
eddy. This process continues till the filament has been 
"packed" as tightly possible. Note that it is the diffusion 
that is responsible for the "packing", i.e., for k > there 
is a limit on how thin a filament can get. Once this sit- 
uation is reached, as the problem is unforced, the whole 
structure remains stationary (or periodic or statistically 
stationary depending on the flow) till the diffusion de- 
stroys all the variance in the scalar field. We argue that 
the eigenmode is a representation of this " packed" struc- 
ture. 



In the eigenvalue formulation, when l s « L ~ l v the 
scalar field is essentially in an infinite domain, the im- 
plication is that the spectrum of C K now typically con- 
sists of a continuous part and hence possesses eigenvalues 
which lie arbitrarily close to each other. As there is no 
dominant eigenvalue we will not see the emergence (at 
short times) of the corresponding eigenmode. 

Taking advantage of l s « l v we linearize Eq. (13) to 
obtain, 



dk(t) 
dt 



-Vu(k ,t) k{t) 



From Eq. (12) we have, 

4>(k ,t) = j>(k ,o) e - K X>»i 2 



(14) 



(15) 



Therefore, in order to see how the moments behave at 
small times we need an estimate of how \k(t)\ evolves. 
Generally this is a fairly intricate question as the solu- 
tion to Eq. (14) involves a time ordered exponentiation 
which reduces to a product of matrices 25 . 



3. The Different Regimes 

Having gotten a feel for the situation when the initial 
scale of variation of the scalar field and that of the 
velocity field are comparable, we proceed to look into 
the evolution of an initially small scale scalar field, i.e., 

l K << l s << l v ~ L. 



1. Useful properties of the FTLE's 

We mention some of the properties of the FTLE's that 
will be needed for calculations in the forthcoming sec- 
tions. In 2D we have the possibility of 2 (asymptotic) 
Lyapunov exponents. Let us denote the larger of these 
by A . Furthermore let us denote the FTLE along a 



particular trajectory, after n iterations, by A(n). Note 
that physically the FTLE's are defined with respect to 
the change in volume of a given set of initial conditions. 
In 2D incompressible systems this only depends on the 
larger eigenvalue in the above mentioned matrix product. 
In higher dimensions the situation is more complicated 
as the FTLE depends on the sum of the positive eigen- 
values (for continuous flows) of the matrix product, i.e, 
it is closely related to the topological entropy of the dy- 
namical system. 

In D dimensions one has the possibility of D (asymp- 
totic) Lyapunov exponents. Let us denote them by 
Aj (i = 1 : D). It is generally assumed that, along a 
particular trajectory, the probability of Aj(n) (i.e., the 
spectrum of Lyapunov exponents calculated after n iter- 
ations) deviating from Aj decays exponentially with n 26 . 
In fact it is generally taken for granted that each A;(n) is 
governed by the central limit theorem (CLT) and is dis- 
tributed around A 27 (see also the discussion in Section 
9.4 of 28 ). We would like to mention that we are not aware 
of a general proof of this statement. Having said this, 
we mention that the results of Balkovsky and Fouxon 7 
appear to lead in this direction though the status of the 
proof of the general statement is not very clear. The 
specific results we are aware of start with the work of 
Furstenberg (see 25 ), who showed that if the matrices are 
i.i.d (i.e. delta correlated in time) random matrices and 
if the system is incompressible, then Ao(= max(A^) > 0. 
A stronger result, which is valid only in 2D, has recently 
been shown by Chertkov et. al. 6 . It states that if the 
tangent maps are random matrices with arbitrary (but fi- 
nite) correlation time then A > 0, moreover A(n) obeys 
the CLT and is distributed around A . 

Therefore, using the CLT for A(n) in 2D, the distri- 
bution of the FTLE's can be expressed as (see Section 
8.6.4 of 29 for the large deviation result), 

Q(A,n) ~ e ~ IlG ( A ~ A o) (16) 

where G(A — Ao) is the Cramer function and Ao is the 
asymptotic Lyapunov exponent. Due to the mixing na- 
ture of the system, as n — > oo the FTLE's along (almost) 
every trajectory will tend to Ao, moreover, the average 
of many realizations along a particular trajectory can be 
interpreted as a spatial average 21 . 



2. The scalar pdf 

At this initial stage the problem is completely reversible 
in the sense that the area of a given blob of scalar is 
invariant. Physically when we release a small blob of 
scalar in a chaotic flow the tendency is for the blob to 
form a filament while respecting the conservation of area. 



Mathematically this implies that on average k(t) grows 
with time, i.e., one of the components of k(t) increases 
while the other decreases. The inference from Eq. (15) 
is that the moments of the scalar field decay in a faster 
than exponential fashion (sec for eg. 30 , 31 ). 

Keeping in mind that the backwards in time problem 
has the same FTLE's as the forward problem, we can 
envision a blob of scalar at a time t to have resulted 
from the diffusive homogenization of a filament (at t = 0) 
which has undergone advective collapse 11 . As the initial 
scale of variation of the scalar field is very small, in effect 
the blob is a result of the diffusive homogenization of 
a large number of independent random concentrations. 
Denoting the scale of the blob by I, for one realization, 
the probability of the blob's concentration being <ft = 4>i 
at time t is 11 , 

P(^) t ~ e -1r s (« (17) 

5(0) is the Cramer function. For notational simplicity 
we have taken the mean of the random concentrations 
making up the blob to be 0. Note that the chaotic nature 
of the flow is embedded in the exponential prefactor to 
S ((/>). For many realizations the average probability of 
the above event is, 

< P(^) t >~ J e - i T^5(0 1 )-tG(A-Ao) dA (lg) 

Here < • > represents an average of many realizations 
along a particular trajectory. As mentioned, this can 
be interpreted as a spatial average. By a steepest de- 
scent argument the value of A that dominates the above 
integral satisfies (denoting it by Ai), 

-fs(0i)e Alt = G'(Ai-Ao) (19) 

As S((t>) > Eq. (19) implies Ai < A . Substituting in 
Eq. (18) we have, 

< P(0!) t >~ e G'(M-Ao) e -Pt . p = G(Ai _ Aq) 

(20) 

Therefore the probability of a particular deviation decays 
exponentially with a rate (3. The important point is that 
(3 is a function of the deviation through Ai. As S(<f>) is 
convex 29 , for <j> = <j) 2 (where \(p 2 \ > \<Pi\), Sfa) > S((j>i). 
From Eq. (19) the dominant FTLE for (f> 2 (denoting it 
by A 2 ) satisfies, A 2 < Ai. Therefore the probability of a 
large deviation is governed by the smaller FTLE's, i.e., 
the tails of the pdf of 4> are sensitive to the tails of the 
FTLE distribution. At the other extreme, as <j> — > 
we have Ai — > A (because S((f>) — > as <f> — ► and 
G'(A - A ) -> as A -> A ). Also, S{<j)) ~ </> 2 for small 
deviations. From Eq. (17) the implication is that, in 
this regime, the scalar pdf has a Gaussian core with a 
width that decreases in time. 



3. Numerical results 

To numerically verify these predictions we use the map, 

x n +i = x n + A sm(By n + ip n ) mod(0, 12n) 
Vn+i =y n + C sin(£>ia; n+ i + a n ) mod(0, 12-k) (21) 

ip-mCn (€ [0, 2ir]) are random numbers chosen at the be- 
ginning of each iteration, B = D\ = i so as to make the 
scale of the flow comparable to the scale of the domain 
and A = 1, C = 1. Diffusion is represented as in Eq. (8) 
with D = 0.25. To satisfy l s << we take the initial 
scalar field to be <fi(x,y) — cos(4x)cos(4y). Fig. 7 shows 
the evolution of the moments for an ensemble average 
over 25 realizations of the map. The faster than expo- 
nential decay of the moments is visible for the first 90 to 
100 iterations. Moreover, the pdf evolution (in steps of 
10 iterations) seen in Fig. 8 goes along with theoretical 
expectations. 



B. Intermediate time scales 

In Fig. 7 there is a transition from faster than expo- 
nential to purely exponential decay of the moments at 
around iteration 100. This transition is a diffusive ef- 
fect (the diffusive time T in the introduction refers to 
the time at which this transition occurs). As mentioned, 
from Eq. (15) it is evident that the growing component 
of k(t) controls the decay of (p(k ai t). In physical space 
the growing component of k(t) represents a shrinking 
physical scale. As there is a lower limit to this shrinking 
scale (i.e. l K ), the magnitude of k(t) saturates when its 
growing component becomes comparable to From 

2 -2 

this time onwards \k(t)\ fluctuates around l K and 
hence, <fi(ko,t) decays in a purely exponential fashion. 
Physically, even though the flow is incompressible, due 
to the nonzero diffusivity the area of a blob of scalar 
is no longer invariant. In fact, the area increases ex- 
ponentially, implying an exponential decay of the scalar 
concentration within the blob itself. 



1. Existing results on a n 

In this regime, one where both l s « L ~ l v and the 
area of a given blob of scalar evolves, the moments 
of the scalar field have been dealt with in detail by 
Balkovsky and Fouxon 7 . Their approach involves shift- 
ing to a comoving reference frame and using the effective 
equation" 



7 6 



The essential point being that now the trajectory equa- 
tion takes the form, 



dt 



= <y aj3 {t) rp 



(23) 



For this problem the Lyapunov exponents arc defined as 
(see for eg. 6 ), 



(24) 



Due to the effective space time separation of the veloc- 
ity field in the comoving reference frame, X(t) is only 
a function of time. In effect the behaviour of A(n) for 
a given trajectory in the FTLE formalism is now valid 
for the entire domain. In fact, the aforementioned re- 
sults of Chertkov et. al. 6 are actually shown in the co- 
moving framework for A(t). Building on this framework, 
Balkovsky and Fouxon 7 derive expressions for the mo- 
ments of the scalar field by considering the evolution of 
an inertia tensor like quantity for a blob of scalar. The 
advantage of the inertia tensor like quantity is that its 
evolution explicitly captures the exponential growth of 
a blob's area which characterizes this regime. The chief 
assumption is that the eigenvalues of the inertia tensor, 
at any finite time t, are governed by the CLT and are 
distributed about the asymptotic Lyapunov exponents 
Xi(i = 1, 2 in 2D). With this they are able to show that 



< \<f>{t)\ n >~ e 



and that a n is a nonlinear function 



of n, implying a non self similar pdf for the scalar field. 



2. A simple example of a single non-self overlapping blob. 

We present a simple example to put some of the notions 
regarding the behaviour of a n in perspective. The ob- 
jective is to analyze the evolution in concentration of a 
single non-self overlapping 8 blob. Let the blob's concen- 
tration be Ct (at time t — T). Now, the increase in area 
of the blob is controlled by the FTLE at the position of 
the blob. Therefore we have, 



t > T 



(25) 



Hence, 



< \C t \ n >= \C T \ n J e-" A (*- T )Q(A,£) dA (26) 



dt 



cr a rf3V a (j) = kV 2 



(22) 



8 What we mean by "non-self overlapping" is a single blob 
that stretches and folds in this large scale velocity field but 
the filaments of the blob do not overlap with each other in 
these intermediate time scales. 



where < • >, as before, indicates an average of many 
realizations over one trajectory. Therefore we have, 

< \C t \ n >= \C T \ n J e -nHt-T)-tG(A-A ) dA ^ 

By a steepest descent argument, for t » T, the above 
integral is dominated for A = A* (implicitly) given by, 



dG, 
dA 



A=A« = ~n 



(28) 



This implies A* < A , moreover for small n Eq. (28) 
implies A* — > Aq. In this region the Cramer function is 
parabolic, i.e., taking G(A — A ) = a(A — A ) 2 we have, 



A* = A - 



n 
2a~ 



Substituting back we obtain (for small n), 



<\C t \ n >^e- Q " t ; a n = n(A 



4a 



(29) 



(30) 



Immediately it is clear that the pdf of the scalar field 
will be non self similar as a n is a nonlinear function of 
n. Furthermore, the nonlinearity is due to fact that there 
is a distribution of FTLE's. If we had a single FTLE or 
if the FTLE distribution collapsed to a delta function in 
a short time then the scaling would be non anomalous 
with a n = uAq. Note that the explicit expression for a n 
given in Eq. (30) is only valid for small n. From Eq. 
(28) the higher moments (as in the previous regime) are 
sensitive to the smaller FTLE's. As our simple example 
follows Balkovsky and Fouxon's work, their form for a n 
(through more detailed considerations which take into 
account the diffusive overlap of many blobs) is similar 
to Eq. (30) (see Eq. 3.9 in 7 ). We reiterate that even 
though the FTLE distribution is responsible for the pdf 
evolution in this regime, fundamentally this is still a dif- 
fusive effect (in the sense that it is the nonzero diffusivity 
that leads to the area evolution which characterizes this 
regime) . 



3. Numerical results 

The upshot of these considerations is that at these in- 
termediate time scales we should observe an evolution in 
the pdf of the scalar field. To verify this, we plot the 
pdf's (in steps of 10 iterations from iteration 90 to iter- 
ation 490) in the upper panel of Fig. 9. The non self 
similarity is clearly evident. The anomalous nature of 
a n can be seen in Fig. 10. We observe that the pdf's are 
characterized by a (progressively smaller) Gaussian core 
and (progressively fatter) stretched exponential tails. In- 
terestingly after about iteration 400 the tails of the pdf 
tend to relax back to a purely exponential form (which 



is the theoretically expected shape, see section IIIB in 
Falkovich et. al. 3 ). The only discrepancy with theoreti- 
cal predictions 7 is that a n does not saturate at high n, 
this might be due to the fact that we are numerically 
unable to get to high enough n to actually observe the 
saturation (note : see the example in the Discussion). 



C. Large time scales 

From the preceeding arguments it would appear that the 
scalar pdf is an ever evolving entity (as is implied in 7 ), 
but once again from Fig. 7 we see that the behaviour 
of the moments undergoes another change between it- 
erations 500 — 600. We claim that it is at this stage 
that the scalar field actually enters the eigenmode, i.e, 
beyond iteration 600 the pdf should become self similar. 
The lower panel of Fig. 9 shows the pdf's (again in steps 
of 10 iterations from iteration 600 to iteration 800). The 
self similar nature of the pdf's is clearly evident, also a n 
(shown in Fig. 10 along with the a n from the previous 
regime) follows a n ~ n. Physically the scale of variation 
of the scalar field is now comparable to that of the veloc- 
ity field, hence the velocity field separation used to derive 
Eq. (22) is invalidated. Correspondingly in the FTLE 
formulation the linearization breaks down. In fact, the 
problem is now in a situation where the mathematical 
considerations leading to the eigenmode (Section II) be- 
come applicable. 



4. Discussion 

In the first part of the paper we described the behaviour 
of a scalar field whose initial scale of variation was com- 
parable to that of the advecting flow (both of which were 
similar to the size of the domain). In both, periodic and 
aperiodic situations we demonstrated the emergence of 
scalar eigenmodes. In spite of the weak barriers that 
exist in the deterministic map used for periodic flows, 
the emergence of a scalar eigenmode is robust (as is seen 
from the two different initial conditions). In the case of 
random maps, representing aperiodic flows, the mixing 
is global (as the barriers are destroyed) and the eigen- 
mode is statistical in nature. Along with the recent work 
of Fereday et. al. 15 this indicates that passive scalar 
advection-diffusion in chaotic flows can be insightfully 
treated as eigenvalue problem as is commonly done for 
steady flows (see for eg. Young et. al. 32 ). 

We then looked at the evolution of a passive scalar 
whose initial scale of variation was small as compared to 
that of the advecting flow (again, the scale of the flow 
was comparable to the size of the domain). Initially, as 



has been noted in previous studies 30 , 31 , the moments of 
the scalar field decay in a faster than exponential fash- 
ion. The pdf of the scalar field in this initial regime is 
an evolving entity, it is shown to be characterized by a 
Gaussian core which shrinks with time. When the scalar 
filaments reach the diffusive scale the behaviour of the 
moments experiences a transition to a purely exponential 
decay. We emphasize that this transition is a diffusive 
effect. In this intermediate regime the pdf of the scalar 
field is still an evolving entity (as is shown by Balkovsky 
and Fouxon 7 , via the nonlinear dependence of a n on 
n). By means of a simple example, the nonlinear nature 
of a n is shown to be dependent on the distribution of 
FTLE's, moreover, the higher moments are seen to be 
sensitive to the tails of the FTLE distribution. Finally, 
when the scalar filamants have stretched and folded to fill 
the domain (or have been "packed" as per our previous 
nomenclature) the field enters an eigenmode. The eigen- 
mode, as before, is characterized by self similar pdf's. 
Numerical results of advection on a lattice followed by 
diffusion appear to confirm these stages of evolution. 

The above mentioned stages that an initially small 
scale scalar field encounters are robust, i.e., once one en- 
ters a particular stage, the statistical properties of the 
scalar field are fixed. A caveat is that the duration of 
the stages, strongly depends on the strength of the ad- 
verting flow. As an illustration, consider the same scalar 
field and map as in Eq. (21), but set A = C = 5 (i.e., a 
flow with stronger stretching properties). The behaviour 
of the moments and a n is shown in Fig. 11 (note that 
in this case, for the intermediate regime, a n does appear 
to saturate for large n) . The stronger nature of the flow 
causes the first stage to be very short, the intermediate 
stage is also relatively shorter as the filaments fill the 
domain quickly. Finally, the eigenmode is realized and 
persists till all the variance is destroyed. 

In essence the picture that emerges is fairly straight- 
forward, as long as there is a valid scale separation (i.e., 
till l s << l v ~ L) the pdf's of the scalar field are evolving 
entities, albeit involving stages characterized by distinct 
decay of moments. Whereas, as soon as l s ~ l v ~ L, 
the scalar field enters an eigenmode with stationary or 
self similar pdf's. Interestingly, some of our prelimi- 
nary numerical work indicates that when we consider 
l v << l s ~ L, the self similarity is maintained whereas 
the shape of the pdf is altered. In fact, smaller the scale 
of the velocity field the more Gaussian is the self sim- 
ilar scalar pdf. As future work we plan to look deeper 
into these cases, especially in light of the recent results 
on non diffusive asymptotic self similarity of turbulent 
decay (i.e. where the advecting flows themselves are self 
similar and nonsmooth) as put forth by Chaves et. al. 33 . 
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Fig. 1. The strong barrier case : The decay of the vari- 
ance for different diffusivities. 





Fig. 2. The weak barrier case. Upper Panels : The two 
initial conditions. Middle Panels : The eigenmodes after 250 
iterations. Lower Panels : The eigenmode after 650 itera- 
tions. 





Fig. 3. The weak barrier case. Upper Panel : The de- 
cay of the various moments (n=2,3,4 and 5 with the higher 
moments appearing lower on the figure) for the first initial 
condition. Middle Panel : The extracted values of a n Vs. n 
(the dash-dot line marked with circles), '+' is a plot of na\. 
Lower Panel : The pdf's in steps of 25 iterations from iter- 
ation 250 to 650 for IC I. The pdf's lie on top of each other 
due to their self-similar nature. 
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Fig. 4. The weak barrier case. Upper Panel : The de- 
cay of the various moments (n=2,3,4 and 5 with the higher 
moments appearing lower on the figure) for the second initial 
condition. Middle Panel : The extracted values of a n Vs. n 
(the dash-dot line marked with circles), '+' is a plot of na.\. 
Lower Panel : The pdf's in steps of 25 iterations from itera- 
tion 250 to 650 for IC II. Once again, the pdf's lie on top of 
each other due to their self-similarity. 
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Fig. 5. The aperiodic (no barrier) case. Upper Panel : 
The decay of the variance. Lower Panel : The higher order 
moments (progressively higher moments appear lower on the 
figure). 




Fig. 6. The aperiodic (i.e., no barrier) case. Upper Panel 
: a n Vs. n (o), the dash dot line is na\. Lower Panel : The 
self similar pdf's. 
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Fig. 7. The decay of the moments for a scalar field whose 
initial scale of variation is small as compared to that of the 
advecting flow (the higher moments are lower on the figure). 
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Fig. 8. The pdf's from iteration 40 to 100 (bold line at 
iteration 100) for the initially small scale scalar field. 
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Fig. 9. Initailly small scale scalar field. Upper panel : 
pdf's from 90 — 490 (non self similar). Lower panel : pdf's 
from 600 - 800 (self similar). 
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Fig. 10. Initially small scale scalar field : a n Vs. n 
(blue (nonlinear) - iteration 90 — 490, red (linear) - iteration 
600 - 800). 




Fig. 11. Case with A=C=5 in Eq. (21). Upper two pan- 
els show the moments (n=2,4,6,8 and 10 with higher moments 
appearing lower on the figure) and the lower panel shows a n 
Vs. n extracted from iterations 20-70 (dashed line) and iter- 
ations 150-300 (solid line) respectively. Note the saturation 
of a n in the intermediate regime. 



